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Abstract 

The vacuum expectation values of the energy-momentum tensor and the fermionic con- 
densate are analyzed for a massive spinor field obeying the MIT bag boundary condition on 
a cylindrical shell in the cosmic string spacetime. Both regions inside and outside the shell 
are considered. By applying to the corresponding mode-sums a variant of the generalized 
Abel-Plana formula, we explicitly extract the parts in the expectation values corresponding 
to the cosmic string geometry without boundaries. In this way the renormalization proce- 
dure is reduced to that for the boundary-free cosmic string spacetime. The parts induced by 
the cylindrical shell are presented in terms of integrals rapidly convergent for points away 
from the boundary. The behavior of the vacuum densities is investigated in various asymp- 
totic regions of the parameters. In the limit of large values of the planar angle deficit, the 
boundary-induced expectation values are exponentially suppressed. As a special case, we 
discuss the fermionic vacuum densities for the cylindrical shell on the background of the 
Minkowski spacetime. 



PACS numbers: 03.70.-Fk, 98.80.Cq, ll.27.+d 



1 Introduction 

It is well known that different types of topological objects may have been formed in the early 
universe after Planck time by the vacuum phase transition [l]. Depending on the topology of 
the vacuum manifold these are domain walls, strings, monopoles and textures. Among them 
the cosmic strings are of special interest. Although the recent observational data on the cosmic 
microwave background radiation have ruled out cosmic strings as the primary source for primor- 
dial density perturbations, they are still candidates for the generation of a number of interesting 
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physical effects sucli as tlie generation of gravitational waves, gamma ray bursts and high-energy 
cosmic rays (see, for instance, [2]). Recently, cosmic strings attract a renewed interest partly 
because a variant of their formation mechanism is proposed in the framework of brane inflation 

In the simplest theoretical model describing the infinite straight cosmic string the spacetime 
is locally flat except on the string where it has a delta shaped Riemann curvature tensor. In quan- 
tum field theory the corresponding non-trivial topology induces non-zero vacuum expectation 
values for physical observables. Explicit calculations for the geometry of a single cosmic string 
have been done for different fields [4]-|29j. Vacuum polarization effects by higher-dimensional 
composite topological defects constituted by a cosmic string and global monopole are investi- 
gated in Refs. [30j for scalar and fermionic fields. Another type of vacuum polarization arises 
when boundaries are present. The imposed boundary conditions on quantum fields alter the 
zero-point fluctuations spectrum and result in additional shifts in the vacuum expectation val- 
ues of physical quantities. This is the well-known Casimir effect (for a review see [31]). In Ref. 
[32], we have studied both types of sources for the polarization of the scalar vacuum, namely, 
a cylindrical boundary and a cosmic string, assuming that the boundary is coaxial with the 
string and that on this surface the scalar field obeys Robin boundary condition. For a massive 
scalar field with an arbitrary curvature coupling parameter we evaluated the Wightman func- 
tion and the vacuum expectation values of the field squared and the energy-momentum tensor. 
The polarization of the electromagnetic vacuum by a conducting cylindrical shell in the cosmic 
string spacetime is investigated in [33] (For a combination of topological and boundary induced 
quantum effects in the gravitational field of a global monopole see Refs. [Ml HH l36].) 

Continuing in this line of investigation, in the present paper we analyze the polarization of the 
fermionic vacuum by a cylindrical shell coaxial with the cosmic string on which the field obeys 
MIT bag boundary condition. We evaluate the fermionic condensate and vacuum expectation 
values of the energy-momentum tensor in both interior and exterior regions of the shell. The 
renormalized vacuum expectation value of the energy-momentum tensor for a fermionic field 
in the geometry of a cosmic string without boundaries is investigated in [TJ |9l [211 [22l [28]. In 
addition to describing the physical structure of the quantum field at a given point, the energy- 
momentum tensor acts as a source of gravity in the Einstein equations and plays an important 
role in modelling a self-consistent dynamics involving the gravitational field. In the problem 
under consideration all calculations can be performed in a closed form and it constitutes an 
example in which the topological and boundary-induced polarizations of the vacuum can be 
separated in different contributions. 

From the point of view of the physics in the region outside the string, the geometry consid- 
ered in the present paper can be viewed as a simplified model for the non-trivial core. This model 
presents a framework in which the influence of the finite core effects on physical processes in the 
vicinity of the cosmic string can be investigated. In particular, it enables to specify conditions 
under which the idealized model with the core of zero thickness can be used. The corresponding 
results may shed light upon features of finite core effects in more realistic models, including 
those used for string-like defects in crystals and superfluid helium. In addition, the problem 
considered here is of interest as an example with combined topological and boundary induced 
quantum effects in which the vacuum characteristics such as energy density and stresses can 
be found in closed analytic form. From the results of the present paper, as a special case, we 
obtain the fermionic Casimir densities for a cylindrical shell with MIT boundary conditions in 
the Minkowski background (for the combined effects of a magnetic fluxon and MIT boundary 
conditions on the vacuum energy of a Dirac field see Refs. [37|). Note that, in addition to tradi- 
tional problems of quantum field theory under the presence of material boundaries, the Casimir 
effect for cylindrical geometries can also be important to the flux tube models of confinement 
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|38[ [32] and for determining the structure of the vacuum state in interacting field theories [ID] . 

We have organized the paper as follows. In the next section, the eigenspinors for the re- 
gion inside the cylindrical boundary are constructed and the eigenvalues of the corresponding 
quantum numbers are specified. These eigenspinors are the basis for the anlysis of the Casimir 
densities in the following sections. In section [3l by using a variant of the generalized Abel-Plana 
formula, we extract from the mode-sum of the fermionic condensate the part corresponding to 
the geometry of a cosmic string without the shell. The part induced by the shell is investigated 
in various asymptotic regions for the parameters. Section [J] is devoted to the investigation of 
the boundary induced parts in the vacuum expectation value of the energy-momentum tensor 
inside the cylindrical shell. The vacuum densities in the region outside the cylindrical shell are 
discussed in section [5l The main results of the paper are summarized in section [HI In appendix 
[3 we give an alternative representation for the fermionic condensate and the expectation value 
of the energy-momentum tensor for a massive fermionic field in the geometry of a cosmic string 
without boundaries. 



2 Eigenspinors inside a cylindrical shell 

We consider the background spacetime corresponding to an infinitely long straight cosmic string 
with the conical line element 

ds^ = dt^ -dr'^ -r'^d4>'^ -dz"^, (2.1) 

where ^ ^ 0o ^ 27r and the spatial points (r, (/>, z) and (r, (/)-|-(/)o, z) are to be identified. The 
planar angle deficit is related to the mass per unit length of the string, //q, by 27r — (/>o = SirGno, 
where G is the Newton gravitational constant. In the discussion below, in addition to the 
parameter ^o, we will use the combination 

q = 27r/0o. (2.2) 

The dynamics of a massive fermionic field is governed by the Dirac equation 

^t'^V^V - mV- = , (2.3) 

with the covariant derivative operator defined as 

= 9^ + . (2.4) 

Here = s'^^-^-y^"'^ are the Dirac matrices in curved spacetime and T^ is the spin connection 
given in terms of the flat space Dirac matrices 7^"^ by the relation 

= l7(«)7We^,)e(,),^^ . (2.5) 

Note that in this formula ; means the standard covariant derivative for vector flelds. In the 
relations above e^^^ is the tetrad basis satisfying ^^g^-^^'ip-^v"'^ — 9^'^ 1 where r/'^'' is the Minkowski 
spacetime metric tensor. 

In this paper we are interested in the change of the vacuum expectation values (VEVs) of 
the fermionic condensate and the energy- momentum tensor for a fermionic fleld, induced by a 
cylindrical shell coaxial with the string on which the fleld obeys MIT bag boundary condition: 

(1 + ^7^^?^^) ^ = , r = a, (2.6) 
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where a is the cyhnder radius and = (0, 1, 0, 0) is the outward-pointing normal to the bound- 
ary. For the evaluation of the VEVs we will use the direct mode summation procedure. In this 
approach we need to have the complete set of the eigenfunctions satisfying boundary condition 

(USD. 

In order to find these eigenfunctions, we will use the standard representation of the flat space 
Dirac matrices: 



7 



(0) 



1 
-1 



, 7 



(a) 



aa 

-(Ta 



, a = 1,2,3, 



(2.7) 



with cTi,(T2,(T3 being the Pauli matrices. We take the tetrad fields in the form used before in 
[H] (see also [281112]): 

/ 1 \ 

fj. _ cos{q(j)) —sm{q4))/r 
W sm{q(p) cos{q(j))/r 

\ 1 / 

where the index a identifies the rows of the matrix. With this choice, the gamma matrices are 
given by 



7° = 7'°', V 



where we have introduced the 2x2 matrices 



e-'"'*' 



-/?' 







~iqcf> 



1 
-1 



For the spin connection and the combination appearing in the Dirac equation we find 



1 



2r 



-7 



and the Dirac equation takes the form 



j'^d^ + -1^1 + ) = 0. 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



For positive frequency solutions, assuming the time-dependence of the eigenfunctions in the 
form e~*'^* and decomposing the bispinor ^ into the upper and lower components, denoted by 
if and Xi respectively, we find the equations 



1 



2r 



P'Ql+^P^ - i ((^ + m) X 



2r 



0, 
0. 



(2.13) 



Substituting the function x from the first equation into the second one, we obtain the second 
order differential equation for the function 99: 



(/p = 0. 



(2.14) 



where n, Z = 1, 2, 3. The same equation is obtained for the function x- 

Because the above equation is in diagonal matrix form, we decompose the spinor 99 into 
the upper, (^1, and lower, 992, components. Taking the eigenfunctions corresponding to these 
components in the form Lpi = Ili(^r)e^^'^^'-'^~^^^~'^^\ I = 1,2, we can see that the solutions of the 



4 



equations for the radial functions, regular on the string, are expressed in terms of the Bessel 
function of the first kind: Ri{r) = ClJp^{\r)^ where the order is defined by the relations 



Pi = \qni + {q- 1)/2|, ^2 = \qn2 - {q - 1)/2|, 
with ni^2 = 0, ±1, ±2, . . ., and 

A = \l iiP' — k'^ — rrfi. 
Hence, the components of the upper spinor are given by the formula 

ipi = CiJi3^ (Ar) exp [i {qnicj) + kz - ujt)] , / = 1, 2. 



(2.15) 
(2.16) 

(2.17) 



Having the upper spinor, we can find the components xi of the lower one by using the first 
equation in (j2.13|) . From this equation we find the following relations 



and 

with the coefficients 



n2 = ni + 1, P2 = Pi + Cni, eni = sgn(ni), 
XI = BiJi3^ (Ar) exp [i {qni4) + kz - Lot)] , / = 1, 2, 



kCi-ienjXC2 „ /cC2 - ie„i ACi 
-Di = ; , = — 



(2.18) 
(2.19) 



u; + m uj + m 

We can see that the bispinor with the components defined by relations ()2.17p and (I2.19P is an 
eigenfunction of the projection of the total momentum along the cosmic string: 



ids + i'T^i^^h^'^^ ] i' = qji^-, 



where 



j=ni + l/2, J = ±1/2, ±3/2, 



(2.20) 



(2.21) 



For the further specification of the eigenfunctions we can impose an additional condition 
relating the constants Ci and C2. As such a condition, following [43| . we will require the 
following relations between the upper and lower components: 



XI = i^Vi, X2 = - — . (2.22) 
From the expressions for the spinor components we find the eigenvalues of the parameter k. 



± sV uP' — k'^ 
k 



, s = ±l. 



and the relation 



C2 = ie„j^(m ± s\/ - k'^)Ci, 
A 



(2.23) 
(2.24) 



for the coefficients in (I2.17p . 

Hence, the positive frequency solutions to the Dirac equation, specified by the set of quantum 
numbers a = {X,j, k, s), has the form 



( ^ft(Ar) \ 

ie,-K,6i+)j^,(Ar)e*^'^ 
K^J/3,(Ar) 



exp ((/(j — l/2)(/) ± fcz — wt)] , 



(2.25) 
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where the orders of the Bessel functions are defined in terms of j as 



Pi = \qj-l/2\=q\j\-ej/2, 
h = \qj + l/2\=q\j\+ej/2. 

In ()2.25p and in the consideration below we use the notations 



A 



(2.26) 



(2.27) 



Note that one has the relation b^g^^ = l/h^ 

The eigenvalues of the radial quantum number A are determined from the boundary condition 
(|2.6p imposed on the eigenspinor (|2.25p . For fixed values of j and s, this leads to the single 
equation 

Jf3^{\a) + ejftW J^2(Aa) = 0. (2.28) 

Using the recurrence relations for the Bessel functions, this equation may also be written in the 
form 

J/3,(Aa) = 0. (2.29) 
Here and in what follows we use the notation 



Jf3^{x) = Xj'p^{x) + (/i - S\/x2 + ^2 _ ejl3i)Jp^{x) 



(2.30) 



with /u = ma. We will denote the solutions of the equation (|2.29p by Aa = A/3^^z, I = 1, 2, . . ., 
assuming that they are arranged in ascending order. Now the set of quantum numbers is specified 
by C7 = {l,j,k,s). 

The coefficient cjj^^ in (|2.25p is determined from the normalization condition 



(2.31) 



where 7 is the determinant of the spatial metric and the integration goes over the region inside 
the cylindrical shell. The delta symbol on the right-hand side of Eq. (j2.3ip is understood as the 
Dirac delta function for continuous quantum numbers {k) and the Kronecker delta for discrete 
ones {l,j,s). Substituting the eigenspinors (|2.25p into Eq. (j2.31|) and using the value of the 
standard integral involving the square of the Bessel function [44J , we find 



Ki + 1 



[2{x^ + /i^) + s (2/3iej + 1) + + 



(2.32) 



with the notation x = Xa. 

The negative frequency eigenspinors are found in a way similar to that used above for the 
positive frequency ones and have the form 



= cl~^ 



exp [-i {q{j + l/2)(j3 + kz - uit)] , 



(2.33) 



iejKsbi'^J^,{Xr)e'i^ 

\ -ie,h^rhp,{Xr)e"i^ J 

with the same notations as in (j2.25p . The boundary condition imposed on this eigenspinor leads 
to the same equation (I2.28P for the eigenvalues of A. The coefficient Ca ^ is found from the 
orthonormalization condition which is similar to ()2.3ip and has the form C, 



(-) 



d^'/bi-\ 
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3 Fermionic condensate 



Fermionic condensate is evaluated by using the mode-sum formula 

(o|vi^|o) = ^V5(-)^(-), 



where ipi- = ipi- is the Dirac adjoint and 



* EE- 

s=±l 1=1 



o- i=±l/2,±3/2,- 

Substituting the eigenspinor ()2.33p into (13. ip . we find 



(3.1) 



(3.2) 



(3.3) 



where Kg is defined by expression (j2.23p and 

P = Pi= q\j\ - ej/2. (3.4) 

The fermionic condensate given by formula (j3.3p is divergent and some regularization procedure 
is necessary. We will assume that a cutoff function is introduced in formula (j3.3p without 
explicitly writing it. 

As the explicit form for X^ i is not known, formula (j3.3p is not convenient for the direct 
evaluation of the condensate. In addition, the separate terms in the mode-sum are highly 
oscillatory for large values of the quantum numbers. A convenient form can be obtained by 
applying to the series over I the summation formula, previously derived in Ref. [35] (see, also, 
|45j ) . In [35] by using the generalized Abel-Plana formula, it has been shown that for a function 
h{z) analytic in the half-plane Rez > and satisfying the condition 



oo 



|/i(2:)| < e(x)e'^l^l , z = x + iy, \z 
with c < 2 and e{x) for x — > oo, the following formula takes place 



(3.5) 



Vr^(A^,;)/i(A^,;) = / hix)dx + -Res 
,1 Jo ^ ^-^ 



vr 



dx 







K^~^^(x) 

-^"/i(xe"/^) f , + e^^*/i(xe-^^/2^^ 



I^-\x) 



(3.6) 



where Ipix), Kfj{x) are the modified Bessel functions, and Tj3{z) is defined by the relation 

v2 



zT^\z) = 4iz) 



sz 



27^2^ 



(3.7) 



In formula (13.61) we used the notations 



zF'jz) + {fi- s^J ii^ -z^ - ejP)F{z) , \z\ < n, 
zF'{z) + (/U =F siy^z"^ - - ej(3)F{z) , \z\ > tJ-, 



(3. 
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for a given function F{z). 

By taking into account the relation 



we can write the mode-sum for the fermionic condensate in the form 

mm = -4;^E^^//5(^' w«)u=A,„ (3.10) 

with the notation 

fpix,y) = (/. - sy^x^ + fi^)jjiy) + (/i + sy^x"^ + ^^)jj^,^{y). (3.11) 
Note that in ([3l^ 

acu = y^x^^^a?T~i^ , (3.12) 

and we have the property 

for X < fi. Now we apply to the sum over / in ()3.10p formula <\3M taking h{x) = xf{x, xr/a)/{auj). 
For this function the residue term in (j3.6p vanishes. The part in the fermionic condensate with 
the last integral on the right-hand side of (13. 6p vanishes in the limit a — > oo, whereas the part 
with the first integral does not depend on a. From here it follows that the latter presents the 
fermionic condensate in the geometry of a cosmic string without boundaries. This can also be 
seen by direct evaluation. 

Indeed, when the cylindrical boundary is absent, the positive and negative frequency eigen- 
spinors are still given by formulae (j2.25p and (j2.33p . where now the spectrum for A is continuous, 
^ A < CO. The corresponding normalization coefficients are found from the condition (|2.3ip 
and have the form 

(CW)-2 = 2vr,^o^(l + &i^)^). (3.14) 



Substituting the eigenspinors into the mode-sum formula ()3.ip . for the fermionic condensate in 
the geometry of a cosmic string without boundaries we find 



A 

to 



j=±l/2,±3/2,-- s=±l 

x[(m - sVA^ + m2)J^(Ar) + (m + s^/ + m'^)jj+,^{Xr)]. (3.15) 

This coincides with the result obtained from the first term on the right of formula (|3.6p applied 
to mode-sum (|3.10p . Formula (j3.15p is further simplified by taking into account the expression 
for (3 and after the summation over s: 



3 

Here and in what follows 



j j=l/2,3/2,- 
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As it is seen from formula (j3.16p . for a massless field the fermionic condensate vanishes in 
the boundary-free cosmic string spacetime. Since the geometry outside of the string is flat, 
the renormalization of the fermionic condensate given by (j3.16|) is done by subtracting the 
corresponding quantity for the boundary-free Minkowski spacetime. The latter is obtained from 
(|3.16p taking in this formula q = 1- Note that for the Minkowski case the summation over j is 
explicitly done by using the formula 

oo 

E['^-i/2(^^) + '^+i/2(Ar)] = 2 Jl{\r) = 1, (3.18) 

j n=0 

where the prime on the sign of summation means that the n = term should be halved. The 
renormalized fermionic condensate in the geometry of a cosmic string without boundaries is 
further investigated in appendix lAl 

From the discussion above it follows that, after the application of the summation formula 
()3.6p . the part in the fermionic condensate with the second integral on the right-hand side 
of this formula corresponds to the VEV induced by the presence of the cylindrical shell. By 
using property ()3.13p and noting that under the change s — > —s we have -F^+-'(x) F^~\x), 
F^~\x) — > F^~^\x), the fermionic condensate in the geometry with the cylindrical shell is 
presented in the decomposed form 

(OIV^^IO) = (0|V;V|0)s + (V5^>cyi. (3.19) 
Here the second term on the right-hand side, 

q \ \ ^ r+°^ r°° dxx 

k''^\x) 

^y— [(/i - is^/^^^^)ll{xT/a) - (^ + is^/^^^^)ll^ {xr/a)], (3.20) 
4 (^) 

is the part in the fermionic condensate induced by the cylindrical boundary. The number of the 
integrations in this formula is reduced by using the relation 

- dxxG{x) TT 

Further, redefining s — > —s in the part of the summation over the negative values j, it can 
be seen that the negative and positive values of j give the same contributions to the fermionic 
condensate. Finally, we arrive to the following formula 



(V'^)cyi = ^-3 ^ / dxxRa 



Ffs {x,xr a) 



(3.22) 



where 

Pj = qj - 1/2, (3.23) 
and the function Fp.{x,y) is given by the expression 



Fp^{x,y) = (^-iVx2-/.2)/|(y) _ (^ + V^2_^2)j|(y). (3.24) 
In (|3.22p . for the modified Bessel functions we use the barred notations 



Fp{x) = xF'p{x) + {fi-Wx^-l^^- P)Fp{x) 

= VFxFp+iix) + Wx^ - fi^)Fpix), (3.25) 
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with F = I,K, and rji = 1, r/x = — 1- 

The boundary induced part (j3.22p is finite for points away from the cyhndrical shell and the 
renormalization is necessary for the boundary-free part, {0\iJjiJj\0)s, only. Note that the ratio in 
the integrand of Eq. ()3.22p may also be presented in the form 



Wfsix) - fi + iy^ 



where 



Wf3{x) = [Ip+i{x)Kp+i{x) - Kp{x)Ip{x)] + 2^1x1 p{x)Kp+i{x). 
For a massless fermionic field, from (j3.22p and (j3.26p we find 



E 



dx X 



I%{xr/a) + ll^_^-^{xr/a) 



(3.26) 



(3.27) 



(3.28) 



In this case the boundary-free part vanishes and the fermionic condensate is always negative. 

Now we turn to the investigation of the fermionic condensate given by Eq. (j3.22p in the 
asymptotic regions of the parameters. For large values of the cylinder radius we use the asymp- 
totic formulae for the modified Bessel functions when the argument is large. By taking into 
account that the main contribution into the integral in ()3.22p comes from the lower limit of the 
integral, to the leading order we find 



(^V')cyl 



^ ^(l + <5(Zj)/J-5/2M 
7 5=±1 



(3.29) 



for am 1. For a massless field, expanding the integrand in (I3.28P we see that the main 
contribution is due to the term with j = 1/2 and one has 



(V'V'>cyl 



2^-iq (r/a)3-i 
'7r2a3 r2((g+l)/2) 7o 



dx 



{xy 



(3.30) 



for r ^ a. 

For points near the string, r/a <^ 1, the main contribution into the boundary- induced VEV 
(j3.22p comes from the lowest mode j = 1/2 and to the leading order we find 



(V'V')cyl 



7r2a3r2((g + l)/2) 



dx x^Re 



'(g-i)/2(^) 

\q-l)/2{x) 



(/i - iy/a 



(3.31) 



This quantity is non-zero in the case of the cylindrical boundary in the Minkowski bulk and 
vanishes for the cosmic string geometry with q > 1. For a massless field formula ()3.3ip is 
reduced to Eq. ()3.30p . For points near the boundary the main contribution comes from large 
values of j and in this case we can use the uniform asymptotic expansions for the modified Bessel 
functions for large values of the order (see, for instance, [IS]). In this way, to the leading order 
we have _^ 

{i^^)cyl ~ - 4^2(a_^)3 ' (3.32) 

for (1 — r/a) <^ 1. As we see, the leading term does not depend on the planar angle deficit and 
corresponds to the same one obtained for a cylindrical boundary in the Minkowski bulk. 

Now we consider the limit when the parameter q is large which corresponds to small values 
of (pQ and, hence, to a large planar angle deficit. In this limit the order of the modified Bessel 
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functions in (|3.22|) is large and we replace these functions by their uniform asymptotic expan- 
sions. Assuming fixed value of the ratio r/a, the integral is estimated by the Laplace method 
and to the leading order we have 

/7/\ ~ g^exp[-(l- (r/a)2)/i/2] (ry 
Hence, for large values of the angle deficit the fermionic condensate is exponentially suppressed. 



4 Energy- momentum tensor 



In this section we consider the VEV of the energy-momentum tensor of the fermionic field inside 
a cylindrical shell on which the field satisfies the boundary condition ()2.6I) . This VEV can be 
evaluated by making use of the mode-sum formula 



(o|r^.|o) = ^ j;[^(-)(x)7(^v,)V^(-)(x) - (V(^V^(-)(x))7.)V'^)(x)] 



(4.1) 



with the negative frequency eigenspinors given by (|2.33|) . We can see that the vacuum energy- 
momentum tensor is diagonal and the separate components are given by the formulae (no sum- 
mation over v) 



(o|r;|o) 



(7 

where the following notations were introduced 



-fp {x,xr/a)\^=xp^v 



(4.2) 



1 



\/x2~+7?y '^^ ' V x"^ + IJ? 



J; 



(y) 



4{y) + 4+eAy) 



2/3 + ej 



J(3{y)Jfi+e, (y), 



(4.3) 



2/3 + ej 

y 



My)Jp+e,{y), ff\x,y) = {kyu')r;'{x,y) 



(0), 



The other notations in (|4.2p are the same as in (j3.1Up . The VEV given by (j4.2|) is divergent and, 
as in the case of the fermionic condensate, it is assumed that a cutoff function is present. Now, 
we can explicitly verify that the VEVs (j4.2|) satisfy the trace relation 

(0|r^^| 0) = m(0|V'V'|0). 

In order to extract from the VEV of the energy-momentum tensor the part corresponding to 
the geometry of a string without boundaries, we apply to the series over / in ()4.2p the summation 
formula (j3.6p with 



h{x) 



— &\x,xr/a). 
auj ^ 



(4.4) 



In a way similar to that used in the case of the fermionic condensate, the VEV can be written 
in the decomposed form: 



r,r 



(4.5) 



where (O |r^| O)^ is the fermionic energy-momentum tensor in the geometry of a cosmic string 
when the cylindrical shell is absent. The second term on the right-hand side of formula (j4.5p is 
induced by the cylindrical shell and is given by the formula (no summation over v) 



cyl 



dx X Re 



IfiAx) 



(4.6) 



11 



where Pj is defined by Eq. (j3.23p . In this formula we have introduced the notations 

= lliy)-ll_,,{y)-{2qj/y)Ip^{y)Ip^+,iy), (4.7) 
F^f{x,y) = {2qj/y)If,^{y)If,^+i{y), 

and Fp'\x,y) = Fp'\x,y). Note that the function for the radial component is also presented 
in the form 

= - (4.8) 

As we see, the vacuum stress along the axis of the string is equal to the energy density. Of 
course, this property is a direct consequence of the bust invariance along this axis. 

By using Eq. (j3.26p . we may write the expressions for the components of the energy- 
momentum tensor in more explicit form given by 

= dxx(i-^v-^) 

Wp^{x)[ll^{xr/a) - /|.+i(xr/a)] - 2/x/|(xr/a) 



for the energy density and by (no summation over v 



Il+iix) + lUx) + 2{fi/x)Ip^{x)Ip^+i{x) 



(4.9) 



IT-) - IX-Tdr 4WpA^)-f,]F^;;\x,xr/a) 
^ -^cyl ^2^42^ ^/|.^^(x)+/|.(x)+2(/./x)/^,,(x)/,^.+i(x)' ^-'^ 



for the radial and azimuthal stresses, ly = 1,2. 

The VEV of the fermionic energy-momentum tensor in the geometry of a cosmic string when 
the cylindrical shell is absent, corresponds to the first term on the right-hand side of summation 
formula (|3.6p . For this VEV we have the mode-sum (no summation over i^) 

(0|T:|0), = -^ E dk dXY, -f^;\\aAr). (4.11) 

The summation over s in this formula is done explicitly and noting that the negative and positive 
values of j give the same contributions, we find 

foo foo \3 



where now 



(0|T;|0), = 4^ / dk dX-g^; {X,Xr), (4.12) 
^ j Jo Jo ^ 



g^f^^i\y) = Jl(y) + jl^^{y)-^j^^{y)jp^^,(y), (4.13) 

4^^^'^) = ^^/3,(y)J,3,+i(y), 5g^(A,y) = (fcV^')4?(^'y)■ 
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As in the case of the fermionic condensate, the VEV (|4.12|) is renormahzed by subtracting the 
corresponding VEV in the Minkowski spacetime without boundaries. The boundary-free renor- 
mahzed VEV of the energy-momentum tensor for a massive fermionic field is further investigated 
in appendix [Xl 

It can be exphcitly verified that the boundary-induced parts satisfy the trace relation {T^)cyi = 
"^(V'V')cyi and the covariant conservation equation Vu{Tp cyi = 0. For the geometry under con- 
sideration the latter is reduced to the single equation 



cyV 



cyl ■ 



(4.14) 



This equation is easily checked by using the relation dy{yF^^\x,y)) = F''^'{x,y) between the 
radial and azimuthal functions in (j4.7p . In the case of a massless fermionic field, for the VEV 
induced by the cylindrical boundary from (|4.9I) and ()4.10p we have (no summation over u) 



(2), 



cyl 



E 



axx — ,n . ^ 6^ 7-^ — [xr/a), 



X 



(4.15) 



where F, 



p ' {y) = Fj^\x,y) for v = 1,2, and the corresponding function for the energy density 
is defined as 



[^{0,0) 



1, 



{y) = -My) 



(4.16) 



In the special case with = 1, from the formulae given above we obtain the fermionic Casimir 
densities for a cylindrical boundary in the Minkowski spacetime. On the left panel of figutre [1] 
we have plotted the corresponding VEVs for the energy density and radial stress as functions 
of the radial coordinate for a massless fermionic field. On the right panel the boundary induced 
parts are presented in the geometry of a cosmic string with the parameter q = 2. 
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Figure 1: Boundary induced parts in the VEVs of the energy density, a^(TQ)cyi (full curves), 
and radial stress, a'^(T/ )cyi (dashed curves), for a massless fermionic field as functions of the 
radial coordinate. The left panel is plotted for q = 1 (Minkowski spacetime) and for the right 
panel q = 2. 

For large values of the cylinder radius the main contribution into the integral in (|4.6|) comes 
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from the lower limit and in the leading order we have (no summation over v) 



_g-2am 



V /cyl 



87ra2 



j 5=±1 



-2am 



, mr ) 



(4.17) 



where = 1,2 and am ^ 1. For a massless fermionic field the main contribution comes from 
j = 1/2 term and from (|4.15p one finds 

2-'?g(r/a)3-i 



for the energy density and 

(^l^>cy, - 



I{q-i)/2{x)K(^g_iy2(.x) - /(g+i)/2(a;)i^{g+i)/2(a;) 







q + 



T 



I V-'O/cyl 



-2q 



cyl 



q + 



T 



l\0 /cyl • 



(4.18) 
(4.19) 



for the radial and azimuthal stresses. 

Now we consider the behavior of the boundary induced VEV of the energy- momentum tensor 
near the string and for points near the boundary. In the first case one has r/a <^ 1 and the 
main contribution in (j4.6p comes from the lowest mode j = 1/2. By using the formulae for the 
modified Bessel functions for small values of the argument, to the leading order one finds (no 
summation over i/) 



q{r/2a) 



q-l 



where 



7rW((g + l)/2)y^ 
/i^/x^ 



fW(. 

F(2)(, 



dx x'?+2Re 



K{q-l)/2{x) 
. I{q-l)/2{x) 



(4.20) 



qF^^\x)=q/{q+l). 



(4.21) 



For a massless fermionic field this formula is reduced to the results given by Eqs. ()4.18p and 

For points near the cylindrical boundary, we replace the modified Bessel functions by the 
corresponding uniform asymptotic expansions for large values of the order. Unlike to the case of 
the fermionic condensate here the leading terms in the VEVs of the energy-momentum tensor 
vanish and it is necessary to take the next terms in the asymptotic expansions. In particular, 
for the function appearing in the integrands of the vacuum stresses we have 



Re 



K^{(3x) 



IpiPx) J 



2f3 7rf 



1-t 



2/3 



(4.22) 



where 



1 



-., r]{x) = \/l + -I- In 



1 + ^ 



(4.23) 



Substituting this into the expression for the azimuthal stress, to the leading order we find 

1 - 5n 



'^yi 607r2a(a 



[l-r/a] < 1. 



(4.24) 
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The corresponding expressions for the radial stress and the energy density are found from the 
conservation equation and the trace relation by using the result (j3.32p for the fermionic conden- 
sate. In this way we obtain the following formulae: 



1 + 10/i 



1207r2a(a 



1 - 5/i 



cyl 



1207r2a2(a_r)2 



(4.25) 



As in the case of the fermionic condensate, the leading terms in the VEVs of the energy- 
momentum tensor components do not depend on the planar angle deficit and are the same as 
for a cylindrical boundary in the Minkowski bulk. Note that, in dependence of the parameter //, 
the vacuum stresses near the boundary can be either positive or negative, whereas the energy 
density remains always negative. As the boundary-free part is finite everywhere outside the 
string axis, for points near the boundary the total VEV is dominated by the boundary induced 
part. 

Similar to the case of the fermionic condensate, it can be seen that for large values of the 
parameter q the boundary induced part in the VEV of the energy-momentum tensor is suppressed 
by the factor (r/a)'^. The dependence of the boundary induced VEVs on the parameter q is 
presented in figure [2] for a massless fermionic field. On the left (right) panel the energy density 
and radial stresses are plotted for the value of the ratio r/a = 0.5 (r/a = 1.5). 
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Figure 2: Boundary induced parts in the VEVs of the energy density, a^(To)cyi (full curves), 
and radial stress, a'^(T/ )cyi (dashed curves), for a massless fermionic field as functions of the 
parameter q. The left panel corresponds to the interior region with r/a = 0.5 and the right 
panel corresponds to the exterior region with r/a = 1.5. 

In figure [3] we give the boundary induced parts in the energy density and the radial stress 
in the geometry of a cosmic string with q = 2 as functions of the mass. The graphs on the left 
panel are for the interior quantities at r/a = 0.5 and the graphs on the right panel are for the 
exterior ones evaluated at r/a = 1.5. As it is seen, we have a non-trivial dependence on the 
mass and in the case of a massive field the polarization effects induced by the boundary can be 
stronger than for a massless field. 
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Figure 3: Boundary induced parts in the VEVs of the energy density, a^(TQ)cyi (full curves), 
and radial stress, a'^{T^) cyi (dashed curves), as functions of the parameter ma for q = 2. The 
left panel corresponds to the interior region with r/a = 0.5 and the right panel corresponds to 
the exterior region with r/a = 1.5. 



5 Vacuum densities in the exterior region 
5.1 Eigenspinors 

In this section we consider the fermionic condensate and the VEV of the energy-momentum 
tensor in the region outside the cylindrical shell. As in the interior case this can be done by 
the direct mode summation. The corresponding eigenspinors have the form (I2.25P and (j2.33p 
with the difference that now, instead of the Bessel functions Jg^ 2('^'^)j linear combinations of 
the functions Jg^ 2('^'^) ^'^d 2(-^'^) should be taken, with Y^{x) being the Neumann function. 
The ratio of the coefficients in this combination is determined from the boundary condition 
(j2.6p imposed on the cylindrical surface. In this way for the positive and negative frequency 
eigenspinors we have the expressions 



( Zfs^{\a,Xr) \ 

iKsenA'^^Zf3^{Xa,Xr)e^i'f' 
KsZi3^{Xa, Xr) 
\ -ienJ^^Zp^{Xa,Xr)e"^'t> J 



exp [±i {q{j =F l/2)(/> + kz - Lot)] 



(5.1) 



where the function Z^{x,y) is defined by the formula 



Zf3^{x,y) = Yp{x)Jf3^{y) - Jp{x)Yp^{y), 



(5.2) 



and 



/?+ = /3i = /?, /3_ = [32, 



(5.3) 



with /3i^2 given by relations p.26p . Here the notation 1^(2;) is defined by p.30p with the 
replacement J — > y. 

The eigenspinors are orthonormalized by condition ()2.3ip . where now the radial integration 
goes over the exterior region. The eigenvalues for A are continuous and on the right-hand side 
of the normalization condition we have 5{X — A'). Since the radial integral diverges for A' = A, 
the main contribution to this integral comes from large values r and we cane replace the Bessel 
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and Neumann functions with the arguments Ar, by the corresponding asymptotic expressions. 
In this way, for the normahzation coefficients in (|5.ip we find: 



1 



(5.4) 



where Kg and bi^^ are defined by Eqs. (j2.23p and (|2.27p . 



5.2 Fermionic condensate 

Substituting the eigenspinors (jS.ip into the mode-sum formula (|3.ip , for the fermionic condensate 
in the region outside the cyhndrical shell we obtain the formula 



(oiv^Vlo) = 



A 



8^'^a;[j2(Aa)+y|(Aa) 
X [(-m + s + m2)Z|(Ao, Ar) - (m + sVA^ + m2)zL^.(Aa, Ar)], (5.5) 



where 



E 



/ + 00 /"OO 



(5.6) 



j=±l/2,±3/2 



s=±l 



In order to extract from the VEV ()5.5p the part induced by the cylindrical shell, we subtract 
the fermionic condensate for the geometry of a string without the shell. As it has been shown 
before, the latter is given by formula (j3.15p . For further evaluation of the difference, we use the 
identities 

^h^rM^^^^) j2 1 V- ■//3(Aa) (02 



Jl{\a)+Y^{\a) 



(5.7) 



(12) 

where cTj = or e^, and Hr '{x) are the Hankel functions. As a result for the fermionic 
condensate we obtain 



(oiv^^lo) 



(0|^V|0)s + 



A Jp{\a) 



x[(m - s^/\^ +m?)Hf'^{\r) + (m + s^/ \^ + m?)Hf^^,{\r)\. (5.8) 



Now, in the complex plane A we rotate the integration contour in the integral over A on the 
right-hand side of formula (j5.8p by the angle tt/2 for / = 1 term and by the angle —tt/2 for I = 2 
term. By using the symmetry properties of the integrands, it can be seen that the parts of the 
integrals over (0, iV /c^ -|- m^) and (0, —i\/k'^ + m?) are cancelled. The number of the remained 
integrations is reduced by using formula (j3.2ip . In this way, introducing the modified Bessel 
functions, we present the fermionic condensate in the decomposed form ()3.19p . where the part 
induced by the cylindrical boundary in the region r > o is given by the expression 



(V'V')cyl 



E 



dx xRe 



Il3i\x] 



---^F'^^f^{x,xr/a) 
(x) 



(5.9) 



where f3j is defined by Eq. (j3.23p and we are using notation 



F'^^f (x, y) = (^ - iV^^^)Kl (y) - (^ + .y^^^)K| (y). 



(5.10) 
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Note that the ratio of the combinations of the modified Bessel functions in (|5.9|) can be written 
in the form 

Ijjjx) ^ W^jx) - n - iy/ x"^ - fJ,"^ 

K^ix) x^[K^^^,ix) + ^ix)]-2fixKpix)Kp+iix)' ^' ' 

where the function Wi3{x) is defined by Eq. (I3.27p . For a massless fermionic field one has the 
formula 

a ^ f°° KjXxr/a) + Kl,^{xr/a) 

and the fermionic condensate is positive. 

Let us consider the behavior of the fermionic condensate in asymptotic regions of the pa- 
rameter. In the limit a ^ with fixed values r, we introduce in (j5.9p a new integration variable 
y = x/a and expand the integrand in powers of a. The main contribution comes from the mode 
j = 1/2 and we have the leading term given below 

^^^^^^^ " vr2r^((g + l)/2)r3 L. - 2"^'^')^(Vi)/2(^) + ^'^(Vi)/2(^)]- (5-13) 

For a massless fermionic field, by using the result from [44J for the integral involving the square 
of the MacDonald function, from here we find 

q{q+ 1) fa\<i 



WW-^^^f^y ' (5-14) 

At large distances from the cylindrical boundary, for a massive field under the condition 
mr ^ 1 the main contribution into the integral in (j5.9p comes from the lower limit of the 
integration and to the leading order we find 



— — / I™ — — 7 r • (5.15) 

VvrH ^ [K j_i/2{ma) 



Here the imaginary part is easily taken by using Eq. (|5.1ip . As we could expect, in this limit 
the VEV is exponentially suppressed. At large distances and for a massless field the behavior of 
the fermionic condensate is described by Eq. (|5.14p . Note that the decreasing of the fermionic 
condensate at large distances is stronger than in the case when the string is absent. For points 
near the boundary, by using the uniform asymptotic expansions for the modified Bessel functions, 
we can see that the leading term in the asymptotic expansion of the fermionic condensate over 
the distance from the boundary is given by the same expression (13.321) as in the interior region. 

5.3 VEV of the energy-momentum tensor 

The VEV for the energy-momentum tensor in the exterior region is found in the way similar to 
that for the fermionic condensate. Here we omit the details of the calculations and give the final 
result. The VEV is decomposed into the sum of boundary-free and boundary-induced parts in 
the form given by (14. 5[) . In the region outside the cylindrical shell the boundary-induced part is 
(no summation over v) 



(5.16) 
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In this formula we introduced the notations 

Fl,f^'\x,y) = Kl{y)-Kl^,{y) + {2qj/y)Kf,^{y)Kp^+,iy), (5.17) 
F^f^'\x,y) = -{2qj/y)Kp^{y)Kp^^,{y), 

and -Fg^^^^^^(x, y) = F^^''^^\x, y). As an additional check we can see that these VEVs satisfy the 
trace relation and the covariant conservation equation. By using formula (j5.1ip . we can write 
the vacuum densities in the form 

Wp^{x)[KUxr/a) - K} ,{xr/a)] + 2^xKI ,{xr/a) 



Kl^.ix) + K^ix) - 2{^,/x)Kp^{x)Kp^+^{x) 
for the energy density and in the form (no summation over v) 



(5.18) 



oo 



x\W(j^{x) — iJi\F^^^'^^''\x,xr /a) 



^^^^^^y' n^a^^J^ "^"^ Kl^,{x)+Kl{x)-2{t,/x)Kp^{x)Kp^^,{xy ^^'^^^ 

for the radial and azimuthal stresses, v = 1,2. We recall that the summation over j in these 
formulae goes in accordance with Eq. (j3.17p . 

In the case of a massless fermionic field from (I5.16P we find the following expressions (no 
summation over v) 

oo 



(r)cvi = J-lTJ dxx'Ft^^'^'\xr/a) 



Ip^{x)Kp^{x) - Ip^+i{x)Kp^+i{x) 
Kl{x) + Kl^,{x) 



(5.20) 



where 



F(;^)(0'0)(y) = h^Kl^^{xr/a) - Kl{xr/a)], (5.21) 

and for the corresponding functions for the radial and azimuthal stresses we haveF^'")(°''^)(y) = 
F(^^)(^)(x,y),z. = l,2. 

Now we turn to the investigation of the VEV in the energy-momentum tensor induced by the 
cylindrical shell in the exterior region in limiting cases. First let us consider the limit a — > for 
fixed values r. Expanding the integrands in powers of a, we can see that the main contribution 
comes from the terms with j = 1/2, and the leading terms are given by the expressions (no 
summation over v) 

, 0\ ^ 2^-'iqm{a/Tf ^ ^ ^ ^ 

^ 7r2r2((g + l)/2)r3 "'''^ [X -mr )K^^_,y,{x), 

{^.)cyl - 7r2r2((c; + l)/2)r3y^/^^ ^(g-l)/2i^'^J' l^-^^i 
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with u = 1,2. For a massless field these terms vanish. In this case from Eq. (|5.2U|) we find the 
following leading behavior (no summation over z^): 



with the coefficients 



An 



cyl 



7r2(q- l)(g + 2)r4 \r. 



q + 3 



-, A2 



-1. 



(5.23) 



(5.24) 



2(g + 4)' g + 4' 

For q = 1 the VEVs behave like {a/r)'^r~'^ln{a/r). 

At large distances from the cylinder and for a massive field the main contribution comes from 
the lower limit of the integral in (j5.16p . By using the asymptotic formulae for the MacDonald 
function for large values of the argument, we find 



cyl 



^(^V)cyl, (T/) 



1 



qme 2™^' \ - 



cyl 



27rr^ 



j 



cyl 



qjRe 



I. 



2mr 



in) 



cyl ' 



ma] 



Kqj-i/2{ma) 



(5.25) 



for mr ^1. As we see, in this limit 



(^/)cyl « (^l)cyl « 



^o')cyr 



At large distances from the 



cylinder and for a massless field the asymptotic behavior of the boundary induced parts is given 
by formula 

The asymptotic behavior of the VEV for the energy-momentum tensor near the cylindrical 
shell is found in a way similar to that for the interior region and the leading terms are given 
by the formulae ()4.24p . (j4.25p . Hence, near the boundary the energy density and the azimuthal 
stress in the interior and exterior regions have opposite signs, whereas the radial stresses have 
the same sign. The boundary induced parts in the VEVs of the energy density and the radial 
stress for the exterior region are plotted in figures [T](3] as functions of the radial coordinate and 
the parameters q and ma. 



6 Conclusion 

In this paper the vacuum polarization effects are investigated for a fermionic field in the geometry 
of a cosmic string with a coaxial cylindrical shell. We have assumed that on the shell the field 
obeys the MIT bag boundary condition. In order to evaluate the fermionic condensate and 
the VEV of the energy-momentum tensor one needs the complete set of normalized eigenspinors 
satisfying the boundary condition. This set for the region inside the cylindrical shell is considered 
in section [2j The corresponding mode-sums for both fermionic condensate and the energy- 
momentum tensor contain series over the zeros of the combination p.30p of the Bessel function 
of the first kind and its derivative. For the summation of these series we used a variant of the 
generalized Abel-Plana formula previously derived in Ref. |35j . This formula allows us to extract 
from the respective VEVs the parts corresponding to the cosmic string geometry without a 
cylindrical shell and to present the part induced by the shell in terms of exponentially convergent 
integrals for points away from the boundary. In this way the renormalization procedure for the 
fermionic condensate and the energy-momentum tensor is reduced to the renormalization of the 
corresponding quantities in the geometry of the boundary-free cosmic string. The renormalized 
VEV of the energy- momentum tensor for a fermionic field in the boundary-free geometry is well 
investigated in literature. In appendix El by using the Abel-Plana summation formula, we give 
alternative integral representations for both fermionic condensate and the energy-momentum 
tensor in the case of a massive field. 
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In the region inside the sheh, the parts in the VEVs induced by the presence of the cyhndrical 
boundary are given by formula (j3.22p for the fermionic condensate and by Eq. (j4.9p . (|4.10p for 
the vacuum energy densities and stresses. These formulae are further simplified for a massless 
fermionic field with the vacuum densities given by Eqs. (j3.28p and ()4.15p . For points near the 
cylindrical shell the boundary induced parts in the VEVs dominate over the boundary-free parts 
and diverge on the cylindrical shell. These type of divergences are well known in quantum field 
theory with boundaries and are investigated for various bulk and boundary geometries. In the 
problem under consideration, the leading terms in the asymptotic expansions in powers of the 
distance from the boundary are given by Eq. (j3.32p for the fermionic condensate and by Eqs. 
()4.24p and (j4.25p for the components of the energy- momentum tensor. These leading terms do 
not depend on the planar angle deficit and are the same as for a cylindrical boundary in the 
Minkowski bulk. The boundary induced parts in the VEVs vanish on the string axis for q > 1 and 
are non-zero in the case of a cylindrical boundary in the Minkowski bulk. Since the boundary- 
free part diverges on the string axis, for points near the string it dominates. For large values of 
the parameter q, which corresponds to a large planar angle deficit, the boundary-induced VEVs 
are suppressed by the factor [r/ay. The boundary induced VEVs have non-trivial dependence 
on the mass of the field and, as it is illustrated by figure [3l for a massive field the polarization 
effects can be stronger than for a massless one. 

Fermionic vacuum densities in the region outside a cylindrical shell with MIT bag boundary 
condition are investigated in section [5l Subtracting from the mode-sums the parts corresponding 
to the geometry of a string without boundaries and by making use of a complex rotation, we 
have derived explicit expressions for the boundary induced VEVs. The corresponding parts 
in the fermionic condensate and the energy-momentum tensor are given by Eqs. (15. 9p . ()5.18p 
and (|5.19p . When the cylinder radius goes to zero, for a fixed value of the radial distance, the 
boundary induced part in the VEV of the energy-momentum tensor vanishes as a'^ for a massive 
field and as a'^^^ for a massless one. At large distances from the cylindrical shell this part is 
exponentially suppressed for a massive field and decay as r~'^{r /aY~^^ in the case of a massless 
field. Note that in the latter case the boundary-free part behaves as and it dominates at large 
distances. For points near the cylindrical shell the leading terms in the asymptotic expansions 
in powers of the distance from the boundary are given by the same formulae as for the interior 
region. In this limit the total VEV is dominated by the boundary induced part. In dependence 
of the mass, the vacuum stresses can be either positive or negative, whereas the energy density 
is positive. In the special case g = 1, from the formulae derived in the present paper we obtain 
the fermionic Casimir densities for a cylindrical boundary in the Minkowski spacetime. 

We have considered the idealized geometry of a cosmic string with zero thickness. A realistic 
model for cosmic string has a non-trivial structure on a length scale defined by the phase 
transition at which it is formed. As it has been shown in Refs. [13 ^ 114 1 [23]. the internal structure 
of the string may have non-negligible effects even at large distances. Here we note that when 
the cylindrical boundary is present, the VEVs of the physical quantities in the exterior region 
are uniquely defined by the boundary conditions and the bulk geometry. This means that if 
we consider a non-trivial core model with finite thickness b < a and with the line element (j2.ip 
in the region r > b, the results in the region outside the cylindrical shell will not be changed. 
As regards to the interior region, the formulae given in this paper are the first stage of the 
evaluation of the VEVs and other effects could be present in a realistic cosmic string. Note that 
from the point of view of the physics in the exterior region the cylindrical surface with MIT bag 
boundary condition can be considered as a simple model of non-trivial string core. 
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A Vacuum densities in the geometry of a cosmic string without 
boundaries 

In this appendix we consider the renormalized fermionic condensate and the VEV of the energy- 
momentum tensor in the cosmic string geometry when the cylindrical shell is absent. For a 
massless field the vacuum energy- momentum tensor was found in [71 E] . Fermionic propagators 
for a massive field are considered in Refs. |21U22j . In the case of a massive field, a representation 
of the VEVs for the energy-momentum tensor in terms of contour integrals is given in [28] . Here 
an alternative integral formulae are given by applying to the corresponding mode-sums the Abel- 
Plana formula. We will do these calculations by using the mode-sum formulae (j3.16p and (j4.12p 
for the boundary-free VEVs. 

First let us consider the fermionic condensate. The renormalization is done by subtracting 
the corresponding quantity for the Minkowski background. The latter is obtained from (|3.16p 
putting q = 1. Substituting in the corresponding formulae 
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integrating over k and A, and introducing a new integration variable y = r"^ /Is^, we find the 
following representation of the renormalized fermionic condensate: 



(V'V'>s,ren = (OlV'^lO)s " (0|V^^|0)m = dye"™ ^ /^^^ 

<5=±1 j 

Next, we apply to the series over j the Abel-Plana summation formula in the form (see, for 
example, [HlllS]) 

E f{n + 1/2) = dx fix) - r dxl^^J^. (A.2) 

It is easily seen that for the summand in (jA.ip the first integral on the right-hand side of (1A.2P 
vanishes and, hence, the divergent parts are explicitly cancelled. Introducing the MacDonald 
function, we arrive to the following expression 

W)ren = ^ / dy ^ '^^'^ \ dx g{q,x)\T^\K,,^y^{y)\. (A.3) 

TT'^'^^ h Jo 

where the notation 

g{q, x) = cosh(7rx) 
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(A.4) 



is introduced. 
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In a similar way we can find the formula for the renormalized VEV of the energy-momentum 
tensor. For the renormalized energy density and the azimuthal stress we have the representations 
given below: 



— 4 ^-OO 

•^0 <5=±1 j 

{T2Xen = V/ dyye-'^-'/'y-yY^6Y^j[q'l^^_s/,{y)-I,.sMy)]. (A.5) 



<5=±1 j 

By making use of summation formula (|A.2p to the series over j, we find the following formulae 
for the renormalized VEVs: 

o poo poo 

(To°)s,ren = -^y^ dy {y + m\^)e-^'''/'y-y dx g{q,x)lm[K,,^y,{y)], 

{Ti)s,ren = dy ye-^'^' '^^y j ^ dxxg{q,x)Re[Ki,+„2{y)]. (A.6) 

The radial stress is found from ()A.3p and (jA.6p by using the trace relation. 

Formulae ()A.6P are further simplified for a massless fermionic field. The integration over y 
is done by using the formula 

^°°....-»K„«/.W = |^P- + 3). (A.7) 
Substituting this into Eq. ()A.6P and integrating over x, we find 

\ U/s,ren g \ ^ /s,ren 28807r r 

In the massless case the radial stress is equal to the energy density. 
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